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Abstract 

We study the application of the classical equations of motion (EOM) within the 
framework of an effective low-energy Lagrangian treated at the loop level. Gauge- 
fixing and ghost terms, which enter naturally in the EOM, are found to lead to 
no physical effects — neither through operator mixing nor in matrix elements. 
Beyond first order in the effective interactions, contact terms have to be included 
when reducing the effective Lagrangian and we present an explicit procedure to 
construct them. Applied to (hadronic) rare B-decays, the EOM drastically sim- 
plify the effective Lagrangian and its matching to the underlying theory, and 
certain cancellations of large (logarithmic) contributions become more transpar- 
ent. Finally, we discuss details of the 'matching' of the effective Lagrangian, 
which may be helpful in incorporating short distance QCD corrections in further 
phenomenological studies. 
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1 Introduction 



Processes in which the external momenta and masses are much smaller than the mass 
scale of the particles that mediate the interactions, are most conveniently described in 
the framework of an effective field theory |T[ . There, the heavy particles are integrated 
out and their relevant effects are summarized in a concise way by local operators in 
an effective Lagrangian. Moreover, short- distance (QCD) corrections can be included 
in a systematic way by means of the renormalization group. Typical phenomenological 
applications of effective theories are quantitative predictions for low-energy processes 
on the basis of a given theory (e.g. the standard model) or when experimental data is 
analyzed in order to obtain bounds or hints for yet unknown 'new physics' (e.g. through 
anomalous gauge-boson couplings). 

The complete set of linearly independent effective operators for certain physical 
processes usually contains many operators which are formally related by the classical 
equations of motion (EOM). It is desirable to exploit these relations among the various 
operators to simplify the expression for a given effective LagrangianQ and at any stage 
of the calculation (e.g. after integrating out heavy particles and before performing the 
renormalization group evolution). To this end, we investigate in this paper the explicit 
form of these relations and their preservation by renormalization effects. 

Of course, at tree (i.e. classical) level, on-shell matrix elements of an effective op- 
erator may be simplified by using the classical EOM, which include the interactions 
from the usual, e.g. QCDxQED, Lagrangian (to account for the fact that internal lines 
entering in the effective operator may actually be off shell). When working at higher 
order in the effective interactions, a more careful formulation of the EOM, which will 
also depend on the effective interactions themselves, is necessary. (Typically, an effective 
Lagrangian is applied only at first order, but one may think of processes, like K° — K° 
or B° — B° mixing, which arise only at second order in the effective weak interactions.) 

At the loop level also renormalization effects have to be taken into account, and they 
are, in fact, crucial for assessing the applicability of the EOM before the renormalization 
group (RG) evolution. The renormalization of composite operators is well-known and 
has been studied in great detail || |], [5|, ||; however, in phenomenological work the issue 
of the EOM is usually either not discussed in a coherent way or the general results are 

1 i.e. a linear combination of effective operators with fixed coefficients. Note the difference to the 
question of simply reducing the set of operators to describe general on-shell amplitudes. For instance, it 
has been stressed || that operators containing covariant derivatives, 1/)^, D^G^ , etc., can be dropped 
from the complete set of operators. However, the explicit form of the relation between the operators 
is not specified; rather their coefficients are to be determined by a complicated matching procedure 
where many additional 1PR graphs must be evaluated. 
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not exploited. 

A complication, which is sometimes overlooked, arises from the fact that the appro- 
priate Lagrangian for the quantized theory contains gauge-fixing terms and Faddeev- 
Popov ghosts. A priori, these terms are to be taken into account in the EOM and they 
can, in principle, contribute to physical processes through mixing with physical opera- 
tors and/or through non- vanishing matrix elements at the loop level. The crucial point 
is that the additional gauge-fixing and ghost terms enter the EOM in a combination, 
which is the BRS variation of some other operator; then, using standard arguments 
]5|, [7|, we find that these additional terms do not affect the matrix elements for phys- 
ical processes. Throughout this paper we shall use the convention that the 'classical 
EOM' include possible gauge-fixing and ghost contributions — in contrast to the 'naive 
classical EOM' without these terms. 

While our discussion of the EOM can be useful in various other phenomenological 
studies (for instance, of anomalous gauge-boson couplings), we focus our attention in 
the second part of this paper on hadronic rare B-decays . In particular, we consider 
here flavour changing neutral b —>■ s currents which arise from loop diagrams, so called 
'penguins', with internal u, c or top quark. The corresponding tree diagrams — if at all 
present — are highly CKM suppressed. 

Penguin diagrams do not only induce interesting exclusive decays, like B — > K*^ 
||, but also contribute considerably to the inclusive hadronic branching ratio for h — > 
s + no charm (via b — > sqq, b — > sgg, etc. M), and they can generate CP violating asym- 



metries in charged B-decays [10]. In various of these processes, one finds a remarkable 
cancellation [T]J of large — since only weakly GIM suppressed — contributions with 
logarithmic dependence on the internal mass. We demonstrate how this cancellation, 
which in fact has been missed in the pioneering works |12[ , is made transparent by 
the EOM. 

The EOM are particularly interesting for calculations of short- distance QCD effects 
(which lead to sizable corrections in the branching ratios [|13| and asymmetries |TJ], [HJ). 
In the abundant literature on such calculations for K- and B-meson decays, essentially 
two approaches are to be found: Many authors perform the RG evolution with an effec- 
tive Lagrangian not reduced by the EOM (see for instance refs. [16, 17, 18]). Of course, 



this is a rather tedious procedure because the corresponding anomalous dimension ma- 
trix is large. In some cases it has been observed that the final result for the process 
considered is independent of whether the EOM are applied before or after performing 



the RG evolution. While in ref. |T7j this fact has been interpreted as a nontrivial result 
of the calculation, the authors of ref. [Q suspected that this is not an accidental fact, 
but did not give a general argument. 
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Other authors, for instance, Grinstein et al. [|T^] and Buras et al. pD| |, employ 
the naive classical EOM to simplify the operator basis right from the beginning. The 
applicability of these EOM is either not discussed explicitly or assumed on the basis 
of a formal proof by Politzer [2l[ (which by a close inspection actually leads to the 
classical EOM including the gauge- fixing and ghost terms). We shall clarify this point 
by demonstrating the appearance of the gauge-fixing and ghost contributions in terms 
of diagrams and by investigating their renormalization. 

Finally, we would like to discuss some details concerning the 'matching' of the full 
theory with the effective Lagrangian (reduced by the EOM), and the recovery of the 
momentum dependent and absorptive parts of the amplitudes in the effective theory 
at leading and next-to-leading log precision. These issues are often somewhat obscured 
in the literature by other technical details related to the actual calculation of the RG 
evolution. 

The remainder of this paper is organized as follows: In section 2 we review the 
proof of the EOM and verify the absence of physical effects from the gauge-fixing and 
ghost terms. We then discuss the appropriate form of the EOM when working at 
higher order in the effective Lagrangian. In section 3 we exploit the EOM to reduce 
the effective penguin operators in b — > s transitions, thereby taking inventory of the 
coefficients with logarithmic mass dependence, and we describe the 'matching' of the 
effective Lagrangian. The discussion in section 4 contains some clarifying remarks, and 
additional details, which may be helpful for further applications, are collected in the 
appendices. 

2 Equations of motion 

2.1 General framework 

The Lagrangian of an effective theory ]l[ consists of two parts: 

C = C Q + £ e ff • 

C Q is renormalizable by power counting and describes those interactions of the particles 
which, in principle, are to be taken into account at arbitrary order in perturbation 
theory. Of course, for a consistent quantization C Q has to include terms for the gauge- 
fixing and Faddeev- Popov ghosts. The effective Lagrangian, C eS , is a linear combination 
of higher dimensional local operators representing the effective interactions. They have 
dimensionful couplings and are to be treated only up to a certain order, usually the first 
one. The effective interactions may be thought as being induced by heavy particles of 
an underlying more fundamental theory. A typical example is £ e g arising from effective 
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weak interactions (induced by virtual W's and top quarks) with C Q = £qcd being the 
usual QCD Lagrangian (neglecting QED interactions for simplicity). 

In order to regularize UV divergences in calculations beyond the tree level we will 
always assume dimensional regularization. The composite operators, Oi(y), in £ e ff are 
defined 0, [7| by their Green's functions (TOi(y)X(xi, . . .)), where T indicates time- 
ordering and the fields are interacting with respect to C D . Physical matrix elements or 
operator identities are obtained by LSZ reduction, and the Green's functions themselves 
are defined perturbatively by (regularized) Feynman graphs. In the following, X will 
always denote a product of elementary fields at different points 

X(xi,...) = ^(x x ) . . . A a ^(x k +i) . . .^(x m ) . . .^f(x n ) . . . . 

For fixed values of the regularization parameter e and of the renormalization scale 
fi, the renormalized effective theory is specified by three elements: (i) The Lagrangian, 
C(g, •••;$), as a function of the fields and parameters (like couplings and masses), 
(ii) the values g R , . . . for these parameters, and (iii) the counter term Lagrangian, C ct 
- or, equivalently, a renormalization prescription to determine C ct from C We call 
£ basic ($) = C(n e g R , ...;$) the 'basic' part of the 'renormalized' Lagrangian £ R (<&) = 
£ basic ($) + C ct (^i e g R , . . . , e; $), where the usual powers of fi e keep the dimensions of 
the g R independent of e. In the following, we will always assume minimal (or ms) 
subtraction for the renormalization prescription. 

Renormalized Green's functions are calculated by using C R and remain finite for 
e — > due to suitable counter terms in C ct . While all counter terms needed for 
Green's functions of elementary fields are proportional to terms in C D (multiplicative 
renormalization of C Q ), the renormalization of insertions of composite operators requires 
additional subtractions proportional to other composite operators. Since C e s is to be 
kept only up to a finite order, only a finite number of such counterterms is needed; 
in this sense C Q + £ e fr is renormalizable. To render the effective operators multiplica- 
tively renormalizable, one includes a priori all interactions in C e R that may appear as 
counterterms. 

The couplings of £ e fj, i.e. the coefficient functions of the operators, are determined 
by 'matching' on-shell amplitudes either with a more fundamental theory or with ex- 
perimental data. By this procedure no distinction is possible between effective theories 
which yield the same on-shell amplitudes. Henceforth, different Lagrangians yielding 
the same on-shell amplitudes will be called 'on-shell equivalent'. An 'on-shell effective 
theory' |I[ can be viewed as a class of on-shell equivalent effective theories. 

If two (effective) theories, renormalized at any scale fi, are on-shell equivalent, then 
the corresponding 'bare' theories^] are also on-shell equivalent, and vice versa. Since 
2 with bare parameters, go(e), . . ., which are divergent functions of e, such that the regularized bare 
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the renormalization group describes the relation between the parameters of the various 
renormalized theories corresponding to the same bare theory but different values of //, 
the RG evolution preserves the on-shell equivalence of effective Lagrangians. Conse- 
quently, one can simplify the expression for £ e s already before the RG evolution by 
replacing C e g by an on-shell equivalent 'reduced' effective Lagrangian. 

To study how the equations of motion can be exploited in order to obtain on-shell 
equivalent effective Lagrangians, we consider a generic effective operator Oeom, which 
contains the fields in a combination that vanishes by the classical EOM derived from 



where $ represents one of the fields from C Q , and Q is a monomial in any of the fields 
and their derivatives (at the same point). Of course, we exclude the case Q = $, since 
then Oeom would correspond to an operator in £ Q and, in particular, would contain the 
kinetic energy term. For Green's functions with an insertion of Oeom one finds || |J 
(after integration over z to avoid derivatives of Dirac delta functions) 



where A[/] = /, d^f, d^d^f, ■ ■ ., denotes the various derivatives (including the fields 
themselves) that may enter in the composite operators Oj. 

For bare Green's functions eq. (§) is readily derived by inspection of the contributing 
Feynman diagrams [see subsection 2.2 for an illustration]. The r.h.s. of (0) consists of 
contact terms which contribute only when z coincides with one of the X\, . . . or yi, . . . , y n . 
They originate from diagrams where the inverse propagator in Oeom cancels either an 
external line from one of the elementary fields in X(xi, . . .), or an internal line ending 
up in a vertex generated by one of the other composite operators 0\(yi) . . -O n {y n ), 
respectively. 

Equation (0) is also valid for the renormalized Green's functions; the derivatives 
5/5$ etc. in (0) and (fj) are then to be taken with respect to the renormalized fields 
and C refers only to the basic part of the Lagrangian. The renormalized version of @ 
is shown either in terms of diagrams M, or using the multiplicative renormalization of 
the composite operators by induction in the number of loops and insertions. 



Green's functions, calculated from the bare Lagrangian £ barc (<i>) = £(g (e), . . . ; <&), differ from the 
renormalized ones only by field-renormalization factors. 



C : 




(1) 




(0 1 (y 1 )...O n (y n )X(x 1 ,...)))dz , 



(2) 
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2.2 EOM at first order in £ e ff 

In many applications of effective theories the interactions from the effective Lagrangian 
are to be treated only at first order. Therefore, we restrict ourselves with the exception 
of subsection 2.5 to this case. We write 

Cf s = g e g c iOi + counterterms , (3) 

where g e s collectively denotes the dimensionful couplings (for instance, g e g = Gf in the 
case of effective weak interactions). All quantities to be considered are at most linear 
in g ef[ or, equivalently, in the coefficient functions q. 

In particular, only Green's functions with single insertions of the composite operators 
from £ e fj are needed, and @ becomes particularly simple 



/(TOeomO) $0i) . . . $(x k )X(x k+1 , . . ))dz 



k 

i . . . ®(xj-x) Q(xj) . . . $(x k )X(x k+1 , . . .)) 

3=1 



(4) 



where $ denotes the field to which the variation in (TI|) refers, and X is a product of 
elementary fields containing no $'s. The contact terms on the r.h.s. do not survive the 
LSZ reduction and, therefore, do not contribute to matrix elements of physical processes 
(see e.g. Joglekar in ref. ]5|, and [|ST| for a discussion of this issue beyond the parton 



level). Hence, any multiple of Oeom can simply be dropped from £ eff . 

Inspecting the reasoning by which @ is derived for the renormalized Green's func- 
tions, one notes that the counterterms for the renormalization of Oeom are themselves 
proportional to operators which vanish by the classical EOM. Of course, the countert- 
erms are in general not proportional to Oeom itself, but in fact related to counterterms 
for Q. For instance, in <3> 3 theory, Oeom — < I )3 (9 2< l ) +m 2 $ + (yf$ 2 /2) requires counterterms 
proportional to <3> 2 (<9 2 $ + m 2 $ + g$ 2 /2) etc. 

The derivation of (0) or (|J) in the path integral representation by a change of vari- 
ables is not very intuitive, and might have even been misleading in cases when ref. |2"T 



is 



quoted in order to justify the naive classical EOM for the gluons (without gauge- fixing 
and ghost terms). To illustrate the derivation in terms of the contributing Feynman 
diagrams we consider here the EOM for the gauge-boson (gluon) field strength 

D^G? = (d,5 ab - g s f abc Al) Or = (J* + Jgf + Jpp): • (5) 

While in the naive classical EOM only the quark contribution 

A a 

{Jqqfa = ~9s 97" y5 > ( 6 ) 

quarks ~ 
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is present, the two other terms on the r.h.s., J g f and Jfp, arise in an unambiguous way 
from the gauge-fixing and ghost terms in the Lagrangian. 

In a covariant gauge (for expressions in a background gauge see appendix A) the 
gauge-fixing and Faddeev-Popov parts of the Lagrangian are given by 

£ gf = -1(&4) 2 , and £ FP = -{dfijDfoh , (?) 

respectively, and the resulting currents on the r.h.s. of (f|) are 

(Jgf)^-V^, and (J FP r a = -g s faUd u r] b )Vc. (8) 

A generic operator vanishing by the classical EOM (0) has the form 

Oeom = Q ■ D^ u - O qq - FP - O gf , (9) 

where Q denotes an arbitrary combination of further fields at the same point and O qq = 
Q ■ Jq q , etc. D^G^ U consists of a one-, a two-, and a three-gluon piece. The one-gluon 
piece, A v — d v d li A tl , can be written in the form Oo + Jgf. Then, J g f cancels against 
the last term on the r.h.s. of @, and the remainder Oo is proportional to the inverse 
gluon propagator {q 2 g^ — ^jrCL^qv in momentum space). 

To illustrate (|j) for Oeom of eq. (|9|), we investigate first the (regularized) diagrams 
which have a vertex from Oq. The external legs of the diagrams are given by the 
elementary fields in @ [with $ = A%] which are on shell for physical matrix elements. 
If the gluon field in Oq is contracted with an external gluon, the diagram contributes to 
the contact terms on the r.h.s. of @. The inverse propagator in Oo acts on the external 
leg and yields zero for on-shell matrix elements. Possible diagrams in which the gluon 
field in Oq is contracted with a gluon in Q at the same point vanish in dimensional 
regularizationP]. In all other diagrams with a vertex from Oo (working at first order in 
£ e ff!), the gluon propagator attached to O ends in an usual QCD-vertex, V, and is 
effectively canceled by the inverse propagator in Oo- 

These diagrams, with a vertex from Oo next to a QCD-vertex V, are in one-to-one 

correspondence to diagrams with one of the remaining vertices from Oeom : A diagram 

with V being a three- or a four-gluon vertex cancels with the corresponding diagram 

containing the two- or the three-gluon piece of D^G^, respectively (see fig. 1). Finally, 

if V is a quark-gluon or a ghost-gluon vertex, then the diagram is canceled by the 

corresponding one with a vertex from O qq or Ofp, respectively. 

3 The vanishing of these 'tadpole'-like diagrams is a crucial property of dimensional regularization 
needed here; otherwise, additional vacuum subtractions are necessary. 
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Similar considerations can be readily carried out for an operator which vanishes by 
the classical EOM for the fermions (quarks), 

ipV =i ^fl - ig^^J ® = m^ , (10) 

or for the case of scalars (see also [0]) and ghosts. 

The main technical complication in proving the renormalized version of (|2|) or (f|) 
by means of diagrams arises from the fact that bare graphs which are in one-to-one 
correspondence (in the sense that they cancel each other) can have a different 1PI 
structure. Thus, their 'forests', which prescribe the renormalization, are not in one- 
to-one correspondence, and it has to be shown that one can remove recursively those 
forests which are not in one-to-one correspondence M. 



2.3 Terms from gauge-fixing and ghosts 

To study the relevance of the terms J g f + Jfp in the EOM for the gauge-boson field 
strength [see the r.h.s. of eq (^)], it is convenient to define two sets of operators: 
(9e : Linear combinations of operators which vanish by the classical EOM derived from 
(the basic part of) £ Q . Of course, terms from gauge-fixing and ghosts in C Q , and the 
EOM for the ghosts themselves [see ( |A.4 )] are to be taken into account as well. 



Ob'- Linear combinations of operators which are the BRS variation of some other oper- 
ator. 

The key observation for the discussion of the gauge-fixing and ghost terms is the fact 
that O g f + Opp belongs to Op provided that Q a v varies under BRS transformations as 

SQ: = gsfa bc Qtv c , (11) 

(with the shorthand notation 5 = Sbrs/Soo, where 5u is the infinitesimal parameter of 
the transformation). Then, 

FP + = 5(-Qld^fj a ) . (12) 



Obviously, eq. ([TT| ) holds if Q^D^G^ is gauge-invariant and contains no ghosts. 
Typical examples for dimension six operators are Q1 = D^G a ^ u or Q a v = ^/ / ~f u X a ^/, 
where \& may be left- or right-handed and a vector in flavour space. After applying the 
EOM twice to {D fl G a (lu ){D\G^ l u ) one arrives at Q a v containing ghosts in the combination 
Qt = (^gf + Jpp)t- m this case an operator from which vanishes by the EOM for 
the ghosts, has to be added on the r.h.s. of QTTD and (TT2]). 

To see that an Opps £ 0~b does not contribute to physical (on shell) matrix elements, 
one writes Oprs = SO and recalls (e.g. 0) 

(T6d(y)X(x 1 , ...)) = -(TO(y)5X{x 1 , ...)). 
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The r.h.s. vanishes after LSZ reduction, because SX contains composite fields which do 
not lead to a physical particle pole. 

It is illustrating to compare the physical matrix elements of O g f, for instance for 
Ql = ^> lu T a ^, with those of J^ a . The fact that matrix elements of [= 5(d"rj a ) G O b ] 
are zero is just the familiar current conservation: The amplitude for u g* — > physical on- 
shell particles" vanishes when one substitutes the polarization vector of the gluon g* by 
its momentum (which can be on or off shell). On the other hand, O g f alone, which is 
not in Ob, has non- vanishing matrix elements. This is readily understood for diagrams 
in which the outgoing gluon line g* from J g f closes a loop (e.g. in fig. 2a): At the 
Ogf-vertex the gluon propagator is indeed contracted by its momentum; however, the 
diagram corresponds to an amplitude where the g* 'decays' into an off-shell particle 
that again enters the O g f-vertex. The overall contribution of these diagrams is canceled 
by the matrix element of Opp in fig. 2b (this can be verified directly by manipulation 



of the corresponding diagrams |22| ). 

As to the renormalization of operators from Ob, one considers an Obrs = 80 G Ob 
and counterterms C, which render all Green's functions of O with elementary fields 
finite. Then, the Green's functions 

(T(0 B rs + SC)(y) ■ X(x u . . .)) R ° = -{T{6 + C)(y) ■ 8X{x x , . . .)) R ° , 

where R indicates counterterms from C , can only be divergent if y coincides with one 
of the Xi (related with a composite field in 5X). Hence, only some operators from Oe, 
together with SC G Ob, are necessary as counterterms for Obrs- 



2.4 Reduction and RG evolution at first order in C e g 

The stability of Oe (respectively Oe © Ob) under renormalization, together with the 
fact that these operators have vanishing on-shell matrix elements, implies that one may 
use the (naive) classical EOM to reduce the basic part of an effective Lagrangian already 
before the RG evolution: If the basic parts of two Lagrangians, renormalized at some 
scale fi, differ only by operators which are in C E (Oe © Ob), then the corresponding 
renormalized theories are on-shell equivalent, and this property is preserved by the RG 
evolution to any other scale. 

To demonstrate this in a more explicit way, we study the renormalization of £ e ff- 
We assume that (0) includes all linearly independent operators which are allowed by 
global symmetries (Lorentz invariance, flavour quantum numbers etc.) and by their 
(canonical) mass dimension. Since Green's functions with single insertions of composite 
operators require only counterterms with at most the dimension of the inserted operator, 
£ e ff can be renormalized multiplicatively by rewriting the effective Lagrangian in terms 
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of renormalized operators^ 

0§ = Z kl OL (13) 

Operator mixing arises when Z is not a diagonal matrix, and we say that an operator 
(or a set of operators) Ok 'mixes into' 0/ if Z ki ^ 0, i.e. if Ok requires counterterms 
proportional to Oi. 

Since C E mixes only into Oe, and since Ob mixes only into Ob © Oe, Z and conse- 
quently the anomalous dimension matrix 

7 = Z± (*->) , (14) 

have block-triangular form, when written in a suitable basis of Oe © Ob and the re- 
maining operators. Therefore, after the RG evolution [] 



.9.0*) f\g) 
exp / ~aFY 9 



CjM , (15) 



ij 



the coefficients of all Oj ^ Oe © Ob do indeed not depend on the initial values 

for the coefficients of the operators from C E © Ob- 

Finally, we note the useful result by Joglekar and Lee || that gauge invariant op- 
erators mix only among themselves and into Oe © Ob- Therefore, the operator basis 
for C e s can be restricted to operators which are linearly independent of Oe © Ob and 
gauge-invariant. Of course, 'unphysical' operators which are not gauge-invariant and 
which are not contained in Oe © Ob can mix into physical operators. However, since 
the latter do not mix into these unphysical operators, their coefficients remain zero if 
the initial values vanish before the RG evolution. This is also illustrated by a recent 



RG calculation of Grinstein and Cho [23[ where Opp was included "for completeness": 
Oep is not in Oe © Ob and indeed mixes into gauge-invariant operators. 



2.5 Equivalent effective Lagrangians at higher order 

When working beyond the linear approximation in the effective interactions one obvi- 
ously has to include the effective interactions themselves in the adequate EOM (and one 
immediately wonders whether to use C Q + C c q or C Q + £ c ff — Oeom in the derivation). 
However, simply removing an operator which vanishes by the classical EOM, does not 
in general lead to a 'reduced'^ effective Lagrangian which is on-shell equivalent to £ e ff- 

4 Alternatively, the counterterms can be viewed as a renormalization of the couplings q; this ap- 
proach is more convenient when working at higher order in £ e g (see appendix B). 
5 The exponential is g-ordered. 

6 Of course, one can always determine a 'reduced' effective Lagrangian £ re d, which is on-shell 
equivalent to £ ff, by performing again the lengthy 'matching' procedure with on-shell amplitudes to 
fix all coefficients for a reduced (e.g. "canonical" [||) set of operators. 
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Green's functions with multiple insertions of Oeom (and £ e ff) have 'non-trivial' contact 
terms which — in contrast to those in (f|) — do not vanish by LSZ reduction. 

In practice, the effective couplings, collectively denoted by g e g, are to be taken into 
account up to a certain order M and the effective Lagrangian has the form C c s = 
g e fi ■ £gg + . . . + g c ff M ■ If £ c fi contains an operator O e "q M , which is of order g^, 

and which vanishes by the classical EOM derived from C D , the 'reduction' of the effective 
Lagrangian by the EOM amounts essentially to the following task: Starting from E q M , 
construct a 'reduced' effective Lagrangian, 

which is on-shell equivalent to C e s up to order g^. 

In the linear case, M — 1, and for m = M one simply has to remove OeSm fr° m 
L e ff. Otherwise, the first step in constructing £ rcd is to extend Oeom recursively to an 
operator 0' EOM which vanishes by the EOM derived form C Q + C' eS . Thereby, C eS is 
defined iteratively by rewriting £ eff as 

C eS = C' cH + 0' EOM . (16) 

The iteration starts with 

r' — r n(n m+l \ 

^eff — ^eff ^EOM ^li/cfT ) > 

o' - o {m) + n (m+1) o(n m+2 ) 

^EOM — ^EOM ^ ^EOM ^ L Ai/cff J ' 

such that 0^q M + O^^P vanishes by the EOM from C Q + C' eS (keeping in C' eS only 
terms of order g eS ) . 

We may assume that E q M contains inverse propagators only through powers of 
5e£ , where 5e = 5/5<& — 8^5/5(3^). Then, 0' EOM can be chosen to contain inverse 
propagators only through powers of 5 E (C Q + C' eS ). For instance, for E q M = Q ■ (5 E C ) 2 
one would set Geom^ = 2Q ■ 5 E C Q ■ S^C^, etc. 

Before removing 0' EOM we have to take into account its on-shell contributions through 
non-trivial contact terms on the r.h.s. of (0) [or its iteration in the case that 0' EOM con- 
tains higher powers of 5e(£ + £g ff )]. The amplitudes are evaluated from the first M 
terms of the Gell-Mann-Low series 

(Texpjz / C eS (y) dy\X(x u . . .)) = ^ — / dy x . . . dy n {TC cS (yi) . . . C eS (y n )X(x u . . .)) 

J n=l n - J 

+ 0{g^)- (17) 

and all on-shell contributions of 0' EOM remaining on the total r.h.s. of flT7|) correspond 
to (sub-) diagrams with the following properties: (a) All internal lines are canceled by 



12 



inverse propagators in the 0' EOM and there is at least one internal line, (b) none of the 
inverse propagators in the 0' EOM ac * s 011 an external line from X(xx, . . .), and (c) all 
vertices from £ Q + C' eS have at least two lines, which are canceled by inverse propagators 
from 0' EOM . 

Obviously, one can just choose suitable subdiagrams with property (a). If (b) were 
not true, the diagram would not survive LSZ reduction; and (c) is a consequence of 
(a) and the assumption that 0' EOM contains inverse propagators only through powers of 

5 E (L + L> eS ). 

Sub-diagrams with property (a) correspond to local vertices which are obtained by 
shrinking the internal lines to a point. In this way one can represent the contribution 
of each contact term by a local operator, and we denote the sum of all these 'contact 
operators' by C[Oe OM ]. More precisely, C[Oeom] = C2 + C3 + ■ ■ ■, where the C n are 
determined by 

n 

i{TC n (y 1 )X(x 1 , . . .))Y[i5(yi - y n ) = 

^t( to eom(2/i) • • • EOM (y n ) exp 

for arbitrary elementary fields, X(x\, . . .), but with the r.h.s. restricted to the contribu- 
tions of maximal connected diagrams satisfying (a)-(c). Working up to order g^, only 
a finite number of contact operators appears in C[Og OM ]. 

Finally, the 'reduced' effective Lagrangian is obtained by replacing 0' EOM by its 
contact operators 

£ re d = C' eS + C[0' EOM } . (19) 

Note that C[Oeom] depends in a non-linear way on 0' EOM (or its coefficient, which we 
have absorbed in O^om)- If £eff = C eS + ^0' EOM , then £ red = C eS + C[XO' EOM ], which 

is given in terms of the contact operators from flTB|) by C[AOeom] = X 2 C<i + X 3 Cs H . 

To illustrate the procedure, we consider a simple example: C Q =—^(d 2 +m 2 )$ 
describes a massive scalar and a massless fermion, and C c s = g c s^^d 2 ^ are 
their (effective) interactions. On-shell amplitudes at order g e g (e.g. one scalar decaying 
into two fermions) are obviously recovered by C^J d = — g eff \I/\I/m 2( I> = £ cff — EOM , 
where EOM = g c s^/^/(d 2 + m 2 )$. At second order in g e s, C e s and ££ d are not on- 
shell equivalent. In particular, C D + C^} d does not yield the correct amplitude for the 
two- fermion scattering, A = —ig 2 g 2 + (s t), where s and t are the Mandelstam 
variables. Following the above procedure, we set 

O'eom = O e 1] om + = 9o^(d 2 + m 2 )$ + g 2 c ^^m 2 , 



and from flT8|), we find C[0' EOM ] = g 2 s /2^^(d 2 + m 2 )(^). One readily verifies that 
£ re d = £ e E — O'eom + ^[^eom] indeed yields the same physical amplitudes as C e ^. In 
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the particular example at hand, this is the case even at arbitrary order in g e g, because 
0' EOM vanishes identically by the EOM derived from C a + C c s — 0' EOM and there are no 
further contact terms at higher orders. 

When working at the loop level, the above reduction procedure remains essentially 
the same: All expressions are to be understood as referring to the basic part of the 
renormalized (effective) Lagrangian, and C[Oe OM ] is still given by the tree result because 
loop diagrams with property (a) vanish in dimensional regularization. Since renormal- 
ized Green's functions of 0' EOM obey the same relations (derived from the basic part 
of the Lagrangian) as the bare ones, is again on-shell equivalent to C^s\ and the 
RG evolution preserves this equivalence. In contrast to the linear approximation in g e g, 
it would be hard to demonstrate this in an explicit way because the reduction by the 
EOM, and the RG evolution (see, for instance, ref. |§) are non-linear operations on C e g. 

For the calculation of amplitudes at higher order in g e $, the Green's functions with 
multiple insertions of C c g are actually needed only after integration over the positions of 
the effective operators [see fllTD ] and not in their general form. Therefore, it is sufficient 
and in fact more convenient to view the counterterms for C c q as a renormalization of 
the effective couplings (see appendix B for more details), rather than of the operators 
as in (|T3|). 



3 Effective treatment of rare B- decays 

We apply now the results of the previous section to rare B-decays; in particular, we 
consider hadronic (and in less detail also radiative) b —>■ s transitions. (The analogous 
b —>■ d modes follow simply from replacing Vt. s by Vtd-) Since the momenta and masses of 
the external particles are at most of the order of the B-mass (and therefore much smaller 
than the mass scales governing the propagation of virtual t-quarks or P^-bosons), it is 
possible to treat these processes in the framework of an effective low energy theory 
The phenomenological procedure in this framework consists of three steps: The deriva- 
tion of the adequate effective Lagrangian, the RG evolution and finally, the evaluation 
of the hadronic matrix elements of the effective operators. 

In step one, the heavy particles are integrated out. The resulting effective action 
contains only the fields of the light particles, however, possibly in a non-local way. In 
order to obtain a local effective Lagrangian one has to perform an operator product 



expansion [p4 in the effective action. If strong (and higher order electroweak) interac- 
tions are absent, this simply corresponds to a Taylor expansion in the external momenta. 
Otherwise, the coefficient functions of the local operators in the effective Lagrangian are 
determined by a 'matching' condition: Appropriate amplitudes calculated in the effec- 
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tive theory are equated with those of the full theory. As soon as divergent loop diagrams 
of effective operators are involved a renormalization scheme has to be specified. The 
coefficients will then depend on this scheme and, in particular, on the renormalization 
scale fi. A natural choice for \x to perform the 'matching' is of the order of the heavy 
masses, My/ in our case. A scale of this size allows to take into account (short-distance) 
QCD effects perturbatively and minimizes logarithmic corrections from higher order 
loops with heavy particles. 

On the other hand, the matrix elements of the effective operators involve logarithms 
of the ratio of [i to the typical mass scale of the low-energy process, which is of order in 
our case. These large logarithms certainly distort a simple perturbative treatment of the 
matrix elements and are removed by the RG evolution in step two: The physical contents 
of the effective Lagrangian, 'matched' at [A ~ Mw, is translated to an effective theory 
with a renormalization scale fi ~ m&. Thereby, the RG evolution of the coefficients [see 
(fjjD] allows to improve the (perturbative) treatment of short- distance QCD corrections 

2 

by summing up all powers of a s £riY^-. 

The RG evolution is straightforward in principle, however, the explicit calculation 
is rather tedious and involves subtle technical details. There are numerous important 
contributions (see, for instance, ref. [2~5| and [17-20]) and we refer to Buras et al. |20 



and Misiak |26[j for the most recent and complete results. Here, we rather discuss 
some details of the derivation and 'matching' of the effective Lagrangian. This provides 
the starting point for further phenomenological studies that intend to include the RG 
improvement. 

The third step, the evaluation of the matrix elements of the effective operators 
with realistic hadron states is, of course, the most difficult task because it requires 
genuine non-perturbative methods (or drastically simplifying assumptions about the 
quark and gluon contents of the hadrons together with some model to describe their 
binding effects). We shall completely ignore non-perturbative effects and stay on the 
parton level throughout. 



3.1 Effective vertices from penguin loops 

In this subsection, we describe one contribution to the effective Lagrangian: Effective 
vertices resulting from one-particle-irreducible (1PI) loop diagrams with internal top 
quark. The application of the EOM already at this stage drastically simplifies further 
calculations (e.g. of the 'full' amplitudes required for the matching). The 'matching' 
procedure to properly define the complete effective Lagrangian is postponed to the next 
subsection; only then, the effective operators corresponding to W-exchange diagrams 
involving only light quarks, like b — > suu and sec, will be introduced (and renormalized) . 
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At leading order in the loop diagrams to be considered are 'penguins' with 

up to three gluons emitted from the t-line inside the loop (diagrams with photons 
and leptons will be discussed below). They are renormalized within the full theory by 
counterterms in the original electroweak Lagrangian (which in fact cancel by the GIM 
mechanism). As long as the penguins do not appear as subdiagrams in two- (or more-) 
loop amplitudes, the external momenta are bound by mb(<C m t or Mw)', and we can 
neglect any momentum dependence which is higher order in m\/{m\ or Mjjp). Then, 
the 1PI diagrams are equivalent to (local) effective vertices and can be summarized in 
terms of a 'penguin' Lagrangian 

, Gp 



-4 



(20) 



We choose the following set of gauge-invariant operators of dimension up to six (see 
appendix C for additional details) 



Pi 

p 3 
Pi 
p 5 
P G 



p 



8(L,R) 



16tt 2 

9s 
16tt 2 

9s 
16tt 2 

9s 
16tt 2 
1 

16vr 2 
1 

16tt 2 
1 

16tt 2 
1 

16tt 2 



s<v (m b R + m s L) f b G£" (= O 

s{iP,a^G^^}Lb, 

s[iP,a^G^^\Lb , 

sip ft ftLb , 

sJDft(m b R + m s L)b , 

si]?M^(L,R)b, 

sMw(m s L,mbR)b , 



(21) 



where 6 and s denote the quark fields, and the covariant derivatives are defined in (||) 
and (jH). 

2 

The coefficients c'j are functions of the variable x = j^- (see appendix C), and c' 7L R 
and c' 8L R depend also on the renormalization scheme of the full theory. Since the b- and 
s-momenta enter P 4 in an antisymmetric way, c' 4 is suppressed by an additional factor 
of the order m^/M^ and can be neglected. We also note that only c' x has a logarithmic 
term, c[ ~ — ^£nx + 0(x), which would dominate for x — > 0. All other coefficients 
approach a constant value in this (unrealistic) limit. 

All coefficient functions are uniquely determined by the (off-shell) digrams for 
the sbg vertex and the sb self energy; a calculation of the 1PI diagrams for b — > sgg 
and b — > sggg is not necessary. In turn, Cp provides the 1PI vertices for all b — > s 
transitions involving up to three gluons (effective operators for more than three gluons 
have dimension higher than six and are suppressed by powers of 1/Myy or 1/mf). This is 
simply a consequence of gauge invariance which is incorporated in £p by using manifestly 
gauge invariant operators. In momentum space the corresponding relations among the 



16 



various effective vertices are much less transparent and the first complete treatments 
of b — > sgg [jO] exploited lengthy Slavnov- Taylor identities or a generalization of Low's 
low energy theorem. 

When Cp enters only at first order, we can take advantage of the EOM (|5|) and 
(|10|), and considerably reduce Cp: P 5 , P 6 and Ptl,r become proportional to the flavour 
off-diagonal mass terms Psl,r with coefficients that cancel when the full theory is renor- 
malized on shellQ Applying the EOM to Pi, one obtains the four-quark operator 

A a A a 
P = s-f v L—b Vl U Y q 

quarks ~~ 

(plus additional operators from gauge-fixing and ghosts, which may be dropped ac- 
cording to the discussion of section 2). Finally, P 3 and P 4 yield color-magnetic dipole 
operators equal to P 2 (P4 yields actually a different chiral structure if m s is not ne- 
glected). 

After this reduction with the help of the EOM, P 2 = O g is the only operator which 
contains a gluon field; and Op and O g are the only two operators remaining in Cp. 
Their coefficients 



9s j 9 2 s 



■(*!(«) + i 



16tt 2 1 167r a ^ iN " y ' 9< 
c' 9 = c 2 + c' 3 = -i(P 2 (x)-| , (22) 



can be expressed in terms of Inami-Lim functions [27 from the sbg- vertex [Pi and P 2 , as 



defined in appendix D, have no constant terms for x — > 0; hence the explicit constants 
in (|2"2"D, which is irrelevant in all physical applications because of the GIM mechanism]. 

Radiative decays require additional penguin diagrams with one, two or three photons 
coupling to the W and to the unphysical Higgs inside the loop. (By choosing a nonlinear 
gauge for the electroweak sector on can avoid diagrams with a photon- py-Higgs coupling. 
Otherwise, these give rise to additional operators ]nj which either vanish on shell, or 
have coefficients that are independent of the internal quark mass and, therefore, cancel 
by to the GIM mechanism.) The resulting effective vertices correspond to four new 
operators, Pi . . . P 4 , which are obtained from Pi ... P 4 by the replacements G a — > F^ u , 
4j- — > 1 and g s — ■> e. Among their coefficients, only c! x has a logarithmic behavior, 
Ci — > — ^£nx, while all others become constant for x — > 0. By applying the EOM to 
Pi ... P 4 (and to P 3 . . . P-j{l,r)i whose covariant derivatives include, of course, also the 
photon field) all local vertices for b — ► #7 or b — > 377 are removed; the only operator 

7 The on-shell renormalization conditions are equivalent to the requirement that all flavour off- 
diagonal mass terms vanish already when the EOM are applied either only to s or only to b. 
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containing a photon field is 7 = P 2 and its coefficient is an other well-known Inami-Lim 
function 

1 ~ 23 

c', = ~c' 2 + ~4 = --(F 2 (x)--) . (23) 

At order a-Gp, additional penguin diagrams with Z-bosons (decaying into leptons), 
and jy-box-diagrams arise. The resulting effective four-Fermi vertices couple sb to 
neutrinos, and to the vector and axial current of the charged leptons. and the 
coefficients of these semi-leptonic operators depend on the gauge that is chosen in the 
electroweak sector. Of course, all (four-Fermi) operators that remain after applying the 



EOM to Pi, have gauge- independent coefficients (see also refs. |27], |28|). 

Cp turns out to be useful even for the treatment of diagrams with light internal 
quarks. Neglecting terms of order -^p- "C 1, where q 2 collectively denotes the external 
momenta, one can always split the amplitude for a penguin diagram with an internal 
quark i = u, c, t into two pieces: 

A ® = E4(#) • (f\Pk\b) tTCC + AA(4) • (24) 

k 1V1 W m i 

For the heavy internal top quark, only the first term, arising from the effective vertices in 
Cp, is relevant. AA vanishes with powers of q 2 /m 2 for momenta small compared to the 
internal quark mass. For light internal quarks the coefficients in the first (local) term are 
evaluated at the appropriate value of m 2 , while the remaining momentum dependence 
of the amplitude is contained in AA (including, for instance, threshold singularities in 
q 2 /m 2 ). 

AA being the difference of the amplitudes at different external momenta is worked 
out most conveniently by using the four-Fermi approximation instead of the full W- 
propagator; no UV divergence arises in A A and no renormalization is necessary. For 
b — > sqq R 

22 

AA(^) = AFx(^) ■ (sqq\Op\b)^ e (25) 

differs from the local part just by a momentum dependent form factor Ai*\ (see ap- 
pendix D). In general, AA consists of matrix elements of local operators, some of which 
may have dimensions higher than six, multiplied by momentum dependent form factors. 



The decomposition (23) is very convenient for processes like b — > sgg etc., where 
the momentum dependence AA contributes significantly in large regions of phase space 
and, hence, can not be neglected. In the effective theory the momentum dependent 
part AA is recovered by loop-level matrix elements of effective four-Fermi operators. 



8 Eqs. (|24j) and ( |25| ) should be contrasted to the more complicated result of ref. |29j| . which keeps 
terms of order q 2 /M^, but is not applicable for large rrii ~ My/. In ( pi| ) the Inami-Lim function from 
c! p gives always the correct mass dependence. 
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However, due to their renormalization, the correspondence is not straightforward and 
has to be clarified in order to incorporate correction factors from short-distance QCD in 
a systematic way. To this end, we take a closer look at the 'matching' conditions that 
define the coefficients of the operators in the effective Lagrangian. 

3.2 Matching of the effective Lagrangian 

In addition to the vertices from £p we consider now the effective operators for the usual 
W-exchange. Writing the effective Lagrangian as 

we need to discuss only C c = Y,k c kOj} in the following; C u is simply obtained from C c by 
replacing the charm field (and mass) everywhere in Ok (and Cfc). In fact, since V u bV* s <C 
V c bV* s , one may neglect C u in most phenomenological applications not concerned with 
CP violation. Note that all terms proportional to VtbV t * (from internal t-quarks) are 
distributed to L c and L u via the unitarity of the CKM matrix. 

The local four-Fermi limit of the usual tree-level VT-exchange, O2 = s^^Lca ■ 
cpj^Lbp, mixes with the operator 0\ = Saj^Lcp ■ cpj^Lbu with reversed color structure. 
Both mix into further four-quark operators, O3 . . . Oq, which carry various color and chi- 
ral structures^, and are related to the penguin operator Op = \[0± + O e — -^(0 3 + O5)]. 
At the two-loop level the four-quark operators mix into the color-magnetic moment 
operator O } 0|. 

According to the discussion of section 2, all operators which vanish by the EOM or 
can be written as BRS variations are irrelevant in the present physical application (first 
order in £ e ffO- m this sense, 0\ . . . 0% and O g form a complete set of the dimension- five 
or -six operators for hadronic b — ► s transitions. In addition, 7 , and the semi-leptonic 
four- Fermi operators, and electromagnetic analogs of O3 ... 0§ have to be included in 
order to describe radiative and semi-leptonic decays, or to take into account mixing at 
order a QE D- 

To obtain a systematic treatment of QCD corrections, the (yet unknown) coefficients 
at fj, — M w are expanded in powers of g 2 s 

Ci {M w ) = cf ] + glc^ + ■ ■ ■ . 

Electromagnetic couplings are included by rewriting e = g s ■ j- and treating j- as an 

independent (running) parameter. 

9 We are here not interested in the explicit expressions for these operators; they can be found in 
e.g. [^0). Note, however, that the Fierz ordering of the operators should be kept fixed throughout all 
calculations in d 7^ 4 dimensions. 
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If the operators are defined with appropriate factors of 1/g 2 p!J they do not mix at 
zeroth order and the expansion of the anomalous dimension matrix starts at order g 2 

— 2 (1) , 4 (2) . 

iki = g s -Tki +9 S -Tki +■■■ ■ 

The RG evolution ([15]) can be written as 

= (4"V) + 9 2 s ■ E$ LL) Qi) + ■■■) c 3 (M w ) , 

where the evolution matrices E sum up all orders of the product a s ■ log{M w / /i 2 ). Since 
the logarithm may be large, this product is not expanded in g s , but rather treated as 
an independent parameter. 

In leading- log (LL) approximation all 0(g 2 ) corrections are neglected, and the of* {Mw) 
are determined by the matching condition 

A «*{b - /) =L cf\M w ) ■ (f\Of\b)^ + {g 2 s ) . (27) 

were the amplitude in the full theory, A iul1 , and the matrix elements, (f\O^\b)^- \ have 
to be evaluated to order g® (or to order g s if A tnU = 0(g s ), like for b — > sg). 
Using appropriate final states, (p7| ) yields 

4 0) = -i , 

4 0) = - U$ ) = -l^(^) + O(^) , (28) 

JO) — fJ I m t \ J ( ml \ _ 1 p / m? \ , ^/ x 
C \ ~ ~ C 7 (^J - -2^(mH + ^(]|r) , 

while all other coefficients are zero at \x = M\y. 

Penguin contributions are order a s in the full theory, and do in LL approximation 
not enter the effective theory at \i = Mw- However, when evolving to /i < Mw, 02 
mixes into O3 . . . 0$ (7^ = — 2^ §§^2 etc.), and if one expands their coefficients to 
leading order in a s , they combine to Op with the coefficient 



ajjj) 2 / fi 2 s (o) 



a s • O(l) + 0(a£) . (29) 



For three-body decays like b — > sgg etc., and for CP- violating asymmetries in charged 
B-decays, one-loop matrix elements of the effective operators are of basic interest be- 



cause they give rise to significant momentum dependent contributions [|11] and to the 
crucial absorptive phase of the amplitudes |L(], [L2|, [J(J . As an example, we consider the 
matrix element 



[sdd\0«\b) 



47T 



srfc/|O P |6) (0) , (30) 
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which corresponds to a penguin diagram with an internal c-quark in the full theory. 
Here, q 2 is the invariant mass of the dd-p&ir and £r = 1 for naive dimensional regular- 
ization with rns subtraction. The one-loop matrix element (^) differs from the pure 
momentum dependence AFi of the full amplitude [see ( p5|) l by additional terms, which 
originate from the renormalization. Eqs. fl29|) and ( p0|) demonstrate how the RG evo- 
lution from jj, = My/ to fi ~ rrib moves the large logarithms from the matrix elements 
into the coefficient functions. 

If the effective theory is matched only with 0(g s ) precision and if the RG evolution 
is performed only in LL approximation, one can not expect that higher order effects, like 
the 0(a s ) matrix element in (|30|), exactly reproduce the full theory. The appearance 
of £r clearly indicates some arbitrariness, because it depends on the details of the 
renormalization scheme. On the other hand, the logarithmic mass dependence and the 
proper momentum dependence is recovered already at this stage: Combining fl29|) and 
(j30l) yields for b — > sdd 



A 



off 



Air 



0(1) 



-in — | 
3 M? 



Q 



m; 



(sdd\O P \b) (0) + 4 o) 



[sdd\O g \b) , (31) 



which is to be compared with the full theory [see (^f)] 



.4 



full 



PCs 

Air 



(sdd\O P \b} i0) + cf ■ (sdd\O g \b) . (32) 



Recalling Fi(x) = — ^inx + 0(x), the two expressions agreeQ within "logarithmic pre- 
cision" ; they differ only by non-logarithmic terms of order a s ■ 1 and, of course, exhibit 
the same momentum dependence. 

To go beyond this precision, one has to include the 0(g 2 ) corrections in the initial 
values of the coefficients and the RG evolution has to be performed in next-to-leading- 
log (NLL) approximation. To this end, the 7^ are needed to evaluate the evolution 
matrix E^ NLL \ and the (Mw) are obtained by solving 

A {ul \b - /) i c! 0) -((/|Of |6)W + g 2 s ■ (f\0*\b) {1) ) + 9 2 s ^-(f\Om {0) + Otf) , (33) 

where A iuU and the matrix elements (f\0*\b) = (f\Of\b)^ + g 2 s ■ (f\Of |&)« + 0{g 4 s ) 
are needed to order g 2 \g\ for A = 0(g s )]. 

For the penguin operators, O3 . . . Og, which arise at fi = My? only in the combination 
Op, the NLL matching condition (0), combined with (^) and (|32"D, yields 



w , 



a s {M ] 



w 1 



47T 



Fii 



m c 



-F X [ 



2 L n ™ c 
+ 3 « E + & M? 



w , 



(34) 



10 ml/M^ and consequently in m^/Myy are considered as O(l); this is consistent with the approx- 
imation that the top and the W are integrated out at the same scale. 
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Similarly, one determines the (scheme dependent) 0(a s ) corrections to the coefficients 
of O12 (see ref. [p0|). Note that the logarithms in (|30| ) and (|34|) combine with AFi{j^j 

and respectively [i = u,c). Therefore, (sqq\Of\b)^ and Cp\Mw) remain 

finite in the limit of negligible internal quark masses (mj — ► 0). 

Since O g and 7 do not mix into the four-quark operators Oy . . . one can con- 
sistently restrict a NLL analysis to this subset [pOfl . At ji ^ Mw, the 0(a s ■ 1) scheme 
dependence entering via the initial values of the coefficients [due to the renormalization 
of the one-loop matrix elements (/ in (|33|)1 is canceled by the corresponding 
scheme dependence in E^ NLL ^ (see Buras et al. [2(| for a detailed analysis of this issue). 

In decay modes, like b —>■ sg, 57, sgg, sg'j, etc., where tree- level and one- loop matrix 
elements contribute at the same (leading) order in g s , all divergencies in the one- loop 
matrix elements of the four-Fermi operators cancel^], and therefore, no renormalization 
scheme dependence remains even in LL approximation. However, a NLL analysis would 
be helpful to reduce the /z-dependenceQ of the physical results, which is the more 
pronounced the larger the QCD corrections are relatively to the uncorrected results at 
fi = M w (e.g. in b — > 37 fTO]). For the operators O g and 7 , the O(g^) matching (P5|) 



'9 ^7' 

requires the finite parts of all corresponding two-loop graphs in the full theory, and the 

(2) 

7y for k = 1. . .6 and I — #,7 have to be determined from tree-loop diagrams. So 
far, such a calculation has not yet been attempted and the mixing of the four-quark 



operators into O g (0 7 ) is carried out only with LL times g s (e) precision [26 



4 Discussion 

In this paper we have studied the validity of the equations of motion within the frame- 
work of an effective field theory described by a Lagrangian £ Q + £ eff , where the inter- 
actions from C (e.g. QCD or QED) are to be treated at arbitrary order and at the 
loop level. We conclude that, in fact, £ e fr can be rewritten by freely using the naive 
classical EOM (i.e. without gauge-fixing and ghost terms) derived from £ Q , provided 
one considers only on-shell matrix elements at first order in £ e fj. In particular, the 
effective Lagrangian may be reduced by the EOM before (RG improved) short- distance 
QCD corrections are evaluated. 

Although this result is not surprising, the reasoning involves non-trivial ingredients, 

11 In fact, this is readily seen with the help of the EOM: The counterterms for e.g. (sgg\02\b)^ are 
proportional to the operator Pi; the latter is equivalent to Op by the EOM and contributes just to the 
renormalization of matrix elements like ( |30| ) . 

12 By construction of the RG evolution, the /i-dependence must cancel order by order in a s . However, 
the RG evolution of the coefficients includes all orders of a s times log, while the matrix elements are 
evaluated only to order g s (g^) in LL (NLL) approximation. 
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which are in principle well-known in a general context [[5], [Tj. The essential steps have 
been carried out explicitly here: Besides renormalization effects for operators which 
vanish by the classical EOM, we have investigated the relevance of the additional terms 
in the correct EOM which arise from the gauge-fixing and ghosts. The possibility of 
dropping these terms is due to the fact that they appear in a combination which is the 
BRS variation of some other operator; such operators do not mix into physical operators 
during the RG evolution and their physical matrix elements vanish identically. 

On the other hand, when the effective interactions from C e g are to be used at higher 
order, operators which vanish by the classical EOM do indeed lead to physical effects. 
We have described an explicit procedure to derive a reduced effective Lagrangian which 
is on-shell equivalent to the original one also beyond the linear order in £ e ff. In this 
procedure an operator Oeom, which vanishes by the classical EOM, is replaced by suit- 
able (physical) operators which account for all non-vanishing contact terms in Green's 
functions with multiple insertions of Oeom- 

Since the renormalization group evolution preserves the relations due to the EOM, 
we note that, it is not necessary for the definition of an ('on-shell') effective theory 
to specify a "canonical form" || of the set of operators: Although the matching with 
on-shell amplitudes determines only a class of on-shell equivalent effective Lagrangians, 
each of them will yield the same physical results after the RG evolution to any other 
scale. 



The ghost operator, Ofp, has also been discussed by Grigjanis et al. |17[| . However, 
there, Opp did not arise through the EOM but was introduced by replacing the 'penguin' 
operator Pi (= O7 in the notation of ref. [|17]]) through a combination of Pi and Opp. 
This ad hoc step was justified by gauge invariance arguments |L7|, Transversality 
is demanded for the sum of the one-loop diagrams that renormalize O7, and Ofp is 
added in order to compensate for a diagram where the two gluons from O7 close a loop 
by emitting one (external) gluon via the three-gluon vertex of QCD. This reasoning 
is not really compelling since only the combination O g f + Opp is BRS invariant (see 
appendix A); in fact, the non-transversal terms in ref. |32] can be compensated by 
counterterms which contain no ghosts but rather vanish by the EOM for the quarks. 
Our discussion shows that Opp arises in a natural way after the use of EOM (|1(]) - 
however, always in combination with O g f. 

The operator O g f was investigated by Eeg and Picek |33[. They claim that the 4% 
piece of the penguin subdiagram gives rise to double-log terms in the mixing of O2 
into 7 , whereas such terms are not reproduced in the effective Lagrangian approach 
TPI , |T7|. We recall that the piece corresponds to the operator O g f and note that 



there exist no two-loop diagrams for b —>■ 37 where Ofp contributes. Thus, our result 
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concerning the absence of physical effects from O g {+Opp implies that the j/qV pieces, and 
the double-logs, in fact cancel for the sum of all diagrams at hand, see fig. 4. (Indeed, 
this can be readily verified by standard Ward-Takahashi identities for the gluon-quark 
vertex. ) 

Besides the striking advantage of the EOM for calculations of short-distance QCD 
corrections, we find the EOM to be particularly useful to keep track of the effective 
(penguin induced) couplings that depend logarithmically on the internal quark masses. 
(Recall that the weak GIM suppression of such terms is, for instance, responsible for 
the amazing fact that a three-body decay like b — > sqq has a larger branching ratio 
than the related two-body decay b — > sg.) In the effective 'penguin' Lagrangian, £p, 
logarithmic coefficients appear only for the operators Pi and Pi, which contain the gluon 
or photon field in the combination D^G^ or d /1 F llu , respectively; via the EOM they 
are equivalent to four- Fermi operators. Consequently, working at lowest order in Gp 
and a s , the amplitudes for any b — > s transition, can receive inx-terms (x = mf/Myy) 
only through contributions from the local four-Fermi operators S'j^L^-b ■ g7 M 4fg or 
s^f^Lb ■ i^i. From this point of view, the crucial — but easily missed — cancellations 



of the inx-terms in b — > sgg etc. JLTJ become rather obvious. 

The distinction between local and momentum dependent parts of the amplitudes, 
which is convenient for purely calculational reasons in the full theory, arises in a natural 
way in the effective theory: There, the two parts originate from tree-level and one- 
loop matrix elements of the effective Lagrangian, respectively, and receive in general 
different corrections when short-distance QCD effects are included. Some processes, 
like b — > sqq, have contributions from 0(1) tree- level matrix elements (of local penguin 
operators) and from 0(a s ) one- loop matrix elements (of usual four- Fermi operators). 
A systematic and unambiguous treatment of both kinds of matrix elements requires 
the use of the next-to-leading results for the QCD-corrections. If only leading-log QCD 
corrections are included, the real parts of the amplitudes can be recovered only up 
to renormalization scheme dependent constants (which, of course, do not affect the 
imaginary parts). Actually, taking the 'leading-log' approximation literally, one can (or 
should) drop the one-loop matrix elements for b — > sqq entirely, and consequently, all 
CP asymmetries become zero [15|]. In other decay modes, like b — > s , ~f,sg,sg'-f, etc., 
no renormalization scheme dependence remains even in LL approximation, but a NLL 
analysis would be desirable to reduce the /x-dependence of the results. 
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Appendices 



Appendix A: Gauge-fixing and ghosts 



In the background field formalism, the gauge field is decomposed into a background field 
B® and the 'quantum field' A a . The Lagrangian derives from the usual Yang-Mills and 
matter Lagrangian by replacing the gauge-field by B° + A a . To fix the gauge, one adds 



£ g f = -2£P?P] , (A.l) 



where we define, for general and $ a , 

Df[V]^ b = d^ a + gf abc V^ c . 
The resulting ghost Lagrangian is 

£fp = fj a D^[B] Dg c [B + A]t] c , (A.2) 
and the corresponding currents on the r.h.s. of @ become 

(Jrf)a = -\ DV aAB] Dl[B] A; , and (J FP )^ = -g s f abc {Dl d [B] fj d ) Vc . 
The theory is invariant under the following BRS transformations 

sa; = D; b [B + A] Vb , 
5b; = o, 

SVa = ~D^[B\Al, (A.3) 
A a 

5# = ig—^Tia , 
S^j = -i g q Va . 

For any Q®, which transforms as 8Q® = gfabcQ^Vc, the operator Q° • (J g f + Jfp)„ can 
be written as a BRS variation 

Q£ • ( Jgf + Jfp); = 5 (-Q» ■ Df[B] fj b ) . 
Since the currents transform as 

$ (JFp)a = gfabc (Jgf + JFp) b Vc , 

5(J g{ r a = -^ ab [B]mB]D c u d [B + A} Vd , 
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the BRS variation of ■ (J g f + Jfp)„ (and of J g f) vanishes upon using the EOM for 
the ghosts 

BZ b [B]I%[B + A]r h = 0. (A.4) 

To obtain the corresponding expressions in a usual covariant gauge, one simply sets 
B^ = everywhere in the above formulae. 



Appendix B: Renormalization at higher order in jC e & 

The non-linear source renormalization (see ref. || for more details) provides a conceptu- 
ally simple method for the renormalization of Green's functions with multiple insertions 
of composite operators. The generating functional 

W = £og J V$ expi J (C Q + J B $ + jfOi)dx , 

includes sources for the elementary fields and for all composite operators (including 
the unity-operator). The bare sources, Jf = Si[jf ,...], are general functions of the 
renormalized sources, Jf(x), and their derivatives. Renormalized (respectively bare) 
Green's functions with insertions of the Oi are given by functional derivatives of W 
with respect to the renormalized (bare) sources. 
The bare sources have an expansion 

$[J* . . .] = (J* + L^d m jfd n jf + 0((J R ) 3 )) [^Z] u , (B.l) 

where the LP^ n are power series of poles in e with residues that depend on the couplings 
of C only; and d m is a condensed notation for possible coordinate derivatives, e.g. 
d° = 1, <9 MlM2 = d 2 /dx^ dx^ , etc. The diagonal matrix // D contains powers of // 
which are needed to leave the mass dimensions of the Jf 1 e-independent and to keep Z 
dimensionless. The matrix of renormalization constants, Z = fi eD Z, is the same as in 



(13|) and provides only the renormalization of Green's functions with a single insertion 
of an Oi. 

To evaluate amplitudes at higher order in £ eff = J^gfOi, onr y Green's functions 
integrated over the coordinates of the insertions of C e s are needed. Therefore, the gf 
can be viewed as additional couplings in the action and their (bare) values follow from 

Since the gf are coordinate independent, only the non-derivative terms (m = n = 0) 
from eq. ( |B.1| ) contribute in Gj. The /^-dependence of the gf is governed by the beta 
functions 





' dG~ 


g? =const. 


[dg R \ 



_1 dG, 
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On the other hand, the RG equation for the (generating functional of) the Green's 
functions 



requires also the knowledge of the derivative terms of eq. QB.1| ) for the 'anomalous 
dimension functions' 

Of course, at first order in the gf, the beta functions of the effective couplings are 
related to the anomalous dimension matrix (14), 



,R\2 



fr^-^gf + OUg 



and the /i-dependence of the "coefficients functions", of eq. ([15]) is equivalent 

to the running of the effective couplings, g~i = q(/i) + 0((g R ) 2 \ governed by the RG 
equation 

[^ + Pw + *wr i(ji ' 9 °' 91 '''' 







Appendix C: Penguin operators 



The operators Pi . . . P$l, as defined in (|21[) , are linearly independent when ignoring 
relations due to the EOM. In fact, together with the analogous operators having opposite 
chirality and with the various four-quark operators, they form a complete set of gauge 
invariant dimension five and six operators for hadronic b — > s transitions. 

The linear independence is readily verified from the Feynman-rules, which we shall 
list here for the b(k) — > s(k) self-energy, written as ^^|£(A;), and for the 1PI vertex 
for b(p) — > s(p') + g(e£,q), which we write as ^^T^p^', q)^-€%- The contributions 
to S and from the various operators are: 



Pi 



Pi 
Pi 

Ps 
Pi 
P 5 
Pe 

PlR 
PsR 









-k 2 fL 
-k 2 (rribR + m s L) 
M^R 
M^m b R 



(p 2 + p' 2 ) 7 M 



{-Tin + 4%) L 
-2ia^ lu q u (m b R + m s L) 

i>Py, - i>p'a ~ p"Pv ~ p"p'n + 2/7^ 



-p 2 + p' 2 ) 7/, + 
(-P 2 7 M " 



11 
P In 



+ PP'n - fan - P"P, 

'hnfi L 

P"ln) ( m b R + m sL) 

M^R 
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The form factors, which enter in a general decomposition of E and T^, 
E(k) = (c M^ + Cl k 2 ) m s L + (d M^ + d^ 2 ) m b R 



r«(p,p',g) 



+ (e M^ + e.k 2 ) J6L + (f Ml + M 2 ) jiR , 

(gooMw + g 01 q 2 + g 02 p 2 + 9o3P 2 ) 7 M 

+ 01 1>Vv + 92 tfp'n + 93 fPv + 94 fp'f, + 95 fit, } 

+ (#67/. i> + 97 i>'lv, + 98Pv + gap'n) • ( m bR + m s L) 
+ h M^ lfl R , 



(C.l) 



(C.2) 



can be found in ref. [|TT|]. They satisfy relations due to Ward-Takahashi identities, 

fx = , (C3) 
and due to the symmetry in p and p' (up to order m\ s /M^) , 



#00 — e o ; h — f , 

C! = di= g 7 + g 9 , a = g i + g 2 + 9i , 

9oi = 1(92- 9x) , 9m = \ {.92 + 9i) + gs 



9o3 = go2 , g3 = g2 , g* = gi , gi = g& , gg = gs ■ 

The remaining form factors are in one-to-one correspondence with the coefficient func- 
tions of the operators P; 



c'i = U92-91) 



208 , 



l(5l+02) 



0, 



4 = 01 + 02 + #5 , 4 = 06 + 08 



(C.4) 



C 7L 



-eo 



-7R 



• 



-8X 



" c j 



-8ii 



-d . 



The relations for photonic operators are similar: The 1PI vertex for b(p) — > s(p') + 
7(6^,9), written in the form -^^T ^(p, p' , q)e fl , vanishes for P x . . . P 4 , while = 
for Pi ... P 4 and = — |T M for P 5 . . . P 8 . The . . . c 4 are obtained by replacing and 
9j in ( |C4| ) by and §j, where the latter are also given in [[IT . 

Appendix D: Inami-Lim functions 

For on-shell quarks (but off-shell gluon) the sbg-vertex has only two independent form 
factors, the Inami-Lim functions P 4 and F 2 



T^(p,p',q) = Fx(x)(q 2 j lx - <£)L - F 2 (x)io ^ u q v (m b R + m s L) . (D.l) 

After applying the EOM to the effective vertices in £p, all coefficients can be expressed 
in terms of Pi and F 2 ; their explicit expressions and the limits for x — > and x — > 1 
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are as follows (77 = -^j) 
FAx) = 



Fo(x) 



J_ 1 3 7_ 2 7 

12 ~ 2^ + 12 77 + 6^" 

2 , 16 « 3 
— te H xtra H — x 

3 3 2 

— -— (x-l) 
72 60 v ; 



1 4 Is 3 '> \ 

77 77 77 mi 

2 ' 3 ' 2 ' ! 



x 



x 



0) 
1) 



-2 V 



-x 



9 



1 



4 



2 

1 1 , 

- + — (x - 1 

8 20 V ; 



4 

(x 

(x ■ 



—77 — -77 + - + -77 + 3i] + -77 £nx 



0) 
1) 



(D.2) 



(D.3) 



For the 567- vertex F\, F 2 have to be replaced by Fi, F 2 in (p.l|) . The magnetic moment, 
F 2 , is given by 

F 2 (x) = ^ 3 + ^rf + ^77 + H + (^-I^ 4 - ^ 3 - ^ 2 - ^77) £nx (D.4) 



7 

— x 
12 

5 13 



x 



0) 



24 120 



(z-1) (a; — ► 1) 



Fx depends on the gauge in the electroweak sector, but combines — via the EOM — with 
the coefficients of the semi-leptonic four- Fermi operators to a gauge-invariant expression 
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The form factor of the momentum dependent part of the amplitude for b — > sqq [see 
]), is given by 



-4/ u(l - u)£n [1 - zu(l - u)\ du 
Jo 



where, setting r = J |1 — A/z 



R{z) 



r+1 



(z<0) 

-2 + z/6 + z 2 /60 + z 3 /420 ■ • ■ (z -> 0) 
2r arctan(r) — m (0 < z < A) 

r in^ 1 + ir tt (z > A) 



(D.5) 
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Figure Captions 



Fig. 1: Subdiagrams with insertions of Oeom of @. The crossed blob <g> denotes O , 
and the hatched ones represent the two- and three-gluon pieces of D^G^ V ' . The dots 
indicate further fields in Q. 

Fig. 2: Diagrams which, in principle, may generate one-loop mixing of O g f + Opp into 
four-quark operators, or to finite contributions to the matrix element for b — > sqq. 

Fig. 3: Diagrams for the mixing of O g f ~ into 7 . The cross x indicates the 
various places where the gluon from O g f has to be attached to. 
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